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Abstract 

We identify the artifact that causes confinement at strong cou- 
pling in compact U(l) lattice gauge theory and show that most of 
the string tension in compact SU(2) lattice gauge theory is due 
to plaquettes of negative trace. 
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1 Introduction 



Virtually all lattice simulations are guided by Wilson's action |]J in which 
the matrices of a compact group play the role of the fields of the continuum 
theory. In these simulations charged particles are confined at strong coupling 
whether the gauge group is abelian or non-abelian |3|, |J . There have been 
a few lattice simulations in which the basic variables are fields. Some of these 
non-compact simulations have no exact gauge symmetry and have shown no 
sign of confinement for either abelian or non-abelian theories [|5], [J . In one of 
them, gauge invariance was partially restored by the imposition of random 
gauge transformations, and a weak confinement signal was observed [0] . But 
it is not clear how much of that confinement signal was due to decorrela- 
tions produced by the noise of the imposed random gauge transformations 
and how much was due to the attractive forces of the gluons. Some very 
interesting simulations || || possess an exact lattice gauge symmetry and 
display confinement for SU(2) and SU(3) theories but not for U(l) above 
(3 = 0.5. The gauge fields of these simulations, however, are not hermitian. 
In all simulations, whether compact or non-compact, confinement has ap- 
peared only when accompanied by significant lattice artifacts. These lattice 
artifacts obscure the confinement mechanising 

We report here the results of two studies of lattice artifacts. The first 
study shows that compact U(l) gauge theory displays confinement at strong 
coupling because its links and plaquettes are circles. One may remove this 
artifact by using a thin barrier to cut the circles of the plaquettes. 

The second study shows that on an 8 4 lattice most of the confinement 
signal in compact SU (2) gauge theory is due to plaquettes of negative trace 
and that, when such plaquettes are avoided, the ratio cx/A| of the string 
tension to the square of the lattice scale parameter is only about 1.6 ± 0.2 
as opposed to about 31 ± 4 when they are included. Since the avoidance of 
plaquettes of negative trace is unlikely to affect perturbative quantities, such 
as the lattice parameter [[□], |i~2l , the decrease in cr/A 2 L may be attributed 
to a decrease in the string tension a. Hence most of the confinement signal 
observable on relatively small lattices in compact lattice gauge theory is due 

§ One may even ask, with Gribov |lC| , whether the lightness of the u and d quarks is 
essential to confinement. 
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to plaquettes of negative trace. 

Since our SU(2) simulations were carried out on a small lattice, we make 
no statement about any possible difference in the continuum limit between 
the physics of the standard Wilson action and that of the positive-plaquette 
Wilson action. In fact we imagine that they are the same. Our result is 
that most of the confinement signal observable with the standard Wilson ac- 
tion on small lattices is due to plaquettes of negative trace, which are lattice 
artifacts. Because the lattice spacing shrinks exponentially as the coupling 
becomes weaker, we expect this result to be valid on lattices as large as 100 4 . 
Quark confinement has only been exhibited on smaller lattices at moderately 
large coupling with actions replete with artifacts. Our result therefore pro- 
vides motivation for further study of confinement on larger lattices at weaker 
coupling with actions that suppress artifacts. 



2 Why Compact U(l) Confines 

Compact lattice simulations of U(l) gauge theory display confinement at 
strong coupling. In Figure 1 we plot the Creutz ratios || we obtained from 
simulations guided by the Wilson action and compare them with the exact 
values of free continuum QED indicated by the curves. At (3 = 0.75, the 
Creutz ratios x(2, 2), x(2,3), x(2, 4), and x(3, 3) all overlap in a striking 
confinement signal. The Wilson action for U(l) confines at strong coupling 
because the links 

U(x, fx) = exp [i9(x, (j,)] = exp [ieaA^x)] (1) 
take values on the unit circle. 

One can avoid this artifact by placing a thin infinite barrier at \8\ = it. 
We used a Metropolis algorithm and rejected any plaquette whose phase 9 
was either greater than 7r — e or less than — 7r + e. We ran on a 12 4 lattice 
and began all runs from a cold start in which all links were unity. Our 
Creutz ratios display no sign of confinement and tend to follow the curves 
of the exact Creutz ratios of free continuum QED. The agreement with the 
exact ratios is better if ones uses an effective inverse coupling (3 e = \ 2 (3 with 
A 2 = 0.73. In Figure 2 we plot the measured and exact Creutz ratios from 
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(3 = 0.25 to /3 = 1.5. Although the continuum theory is free and exactly 
soluble, the lattice action is non-linear and does require renormalization. 
The renormalization g e = go/X is the simplest rule that in the continuum 
limit satisfies 

Jim \9 2 ~9o\ = Q (2) 

90-0 g Q 

We saw no confinement signal as long as the step size was smaller than the 
thickness 2e of the cut in the phase 9. We took .024 < e < 0.100 and noticed 
no sensitivity to e within that range. The nearly overlapping points at (3 = 1.5 
correspond to e = .024 and to e = 0.1. We also performed simulations with 
Manton's action using e = 0.026 and found a similar absence of confinement. 

Confinement signals arise when links are decorrelated. At strong coupling 
many U(l) plaquettes are near —1. Thus when a link on a plaquette that is 
near —1 is being updated, that link can jump to a value (often far from 1) 
that pushes the plaquette past the point —1 of maximum action to a lower 
action. After many such events, the links are decorrelated, and a confinement 
signal appears. 

One may also interpret these results in terms of monopoles HT5| . When 
the phase 9 of each plaquette is required to lie between 7r — e and — n + e, it 
follows that for e > 0, no string can ever penetrate any plaquette. Although 
the circle artifact is frequently so interpreted, our purpose here is rather to 
show that by removing it, we may improve the agreement between continuum 
and lattice QED. 



3 Compact 5/7(2) 

In view of these results for U(l), one might wonder whether similar lattice 
artifacts exist in the case of the group SU(2). Because U(l) and 577(2) have 
different first homotopy groups (7i"i(Z7(l)) = Z but 7i"i (577(2)) = 0), it would 
seem likely that the excision of a small cap around the antipode (g = — 1) 
on the 577(2) group manifold (the three sphere 53 in four dimensions) would 
have little effect upon the Creutz ratios. 

To check this assumption, we ran from cold starts on an 8 4 lattice and 
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used a Metropolis algorithm with a modified form of the Wilson action in 
which we rejected plaquettes that lay within a small cap around the antipode. 
Specifically the allowed plaquettes P had to satisfy the rule 

]f£iP > - cos(^r) « -0.951. (3) 
In terms of the parameterization 



exp I i— ■ & I = cos — / + id ■ a sin — , (4) 

I 2 J I 2 / \ 2 ) 

the excluded cap is described by the condition 

— <\e\< —. (5) 

5 - 1 1 - 5 V ; 

The step size was small compared to the size of the excluded cap. As ex- 
pected, the Creutz ratios x(hj) do not depend upon whether the small cap 
was excluded. At (3 = 2, for example, the Creutz ratios with the cap in- 
cluded are: x(2,2) = 0.5995(1), x(2, 3) = 0.5813(4), x(2,4) = 0.582(5), 
x(3,3) = 0.556(5), and x(3,4) = 0.54(2); while with the cap excluded they 
are: X (2,2) = 0.596(4), X (2,3) = 0.574(7), X (2,4) = 0.57(1), X (3, 3) = 
0.47(5), and x(3,4) = 0.5(1). The exclusion of the cap makes little differ- 
ence. 

But what about the excision of a large cap? To study this question, we 
performed simulations using a Metropolis algorithm and a modified form of 
the Wilson action in which we rejected all plaquettes that had a negative 
trace, thus excluding half of the SU{2) sphere. This positive-plaquette Wil- 
son action has been studied by other groups (JT2L TA\. We used an 8 4 lattice 



and began with all SU(2) group elements equal to the identity (cold starts). 
We used a very small step size; the maximum change in any of the four real 
numbers that describe each group element was 0.005. Because of the small 
step size, we allowed 2,000,000 sweeps for thermalization. 

In our positive-plaquette simulations, the Creutz ratios x(hj) as func- 
tions of (3 tend to follow the values given by the SU{2) tree-level perturbative 
formula [|T5| 

3 

Xo{i,j,P) = [-u(i,j) -u(i - 1) +«(«', j - l)+u(i - (6) 
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where 



u(i, j) = - arctan - + - arctan - — log I — H — - I . (7) 

J 3 1 1 




The fit of the perturbative formula to the data improves if we use an 
effective inverse coupling (3 e = \ 2 f3 which corresponds to a renormalized 
coupling constant g = go/X with < A < 1. This scheme offers the simplest 
renormalization of the coupling constant that in the continuum limit satisfies 
the condition P). In Figure 3 we plot the Creutz ratios x(hj) we obtained 
in simulations guided by the positive-plaquette Wilson action. The curves 
display the perturbative values as given by the tree-level formulas @ and 
(0) with the effective inverse coupling (3 e = 0.54 (3. The data follow the 
perturbative curves near (3 = 3 but not near (3 = 1 where the bunching of 
the Creutz ratios for different i and j indicates confinement. 

In the confining phase, the Creutz ratios x(hj) ought to be given approx- 
imately by the product of the string tension a and the square of the lattice 
spacing a 2 ((3 e ). For the group SU(2), the two- loop result for the dependence 
of this product aa 2 ((3 e ) upon the effective inverse coupling (3 e is 



aa 2 ((3 e ) ta exp 



6n 2 p e 102, [6ir 2 p e 



11 121 ° 11 



(8) 



where is the lattice scale parameter 0. We expect this scaling formula 
to hold in a transition region where perturbation theory is still valid and 
where the quark- ant iquark static potential V(r) is a linear combination of a 
confining potential and a Coulomb potential: 

ft = Tr T [^ max ~ ^) Cr ° 2 (/ 3 e) + (P~ (3min)Xo(i,j,f3e) (9) 

in which the string tension aa 2 ((3 e ) is given by the scaling formula (§) and 
the perturbative Creutz ratio Xo(hj,Pe) is given by the tree-level formula 
(|^ 0). Near one end of this region, the potential V(r) is mostly Coulomb; 
near the other end, V(r) is mostly linear. 

In simulations guided by the positive-plaquette Wilson action, we found 
a wide transition region between (3 = 1.0 and (3 = 3.0. Our best fit to the 
Creutz ratios X (2,2,/3), x(2,3,/3), x(2,4,/3), and x(3,3,/3) at (3 = 1.0, 1.25, 
1.5, 1.75, 2.0, 2.25, 2.5, 2.75, and 3.0 was (3 e = 0.36/5, (3 min = -2.77, /3 max = 
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6.51, and appj A| = 1.61±0.2. In the interpolation ([|), the coefficient (/? max — 
/3)/(/3 max — /3 min ) of the positive-plaquette string-tension term appa 2 (/3 e ) runs 
from 0.59 at (3 — 1.0 to 0.38 at (3 = 3.0, and the coefficient (/3 — /3 m i n ) /(/3 max — 
Anin) of the tree-level perturbative term Xo(hj, A) runs from 0.41 at j3 — 1.0 
to 0.62 at (3 — 3.0. Our fit to the data is displayed in Figure 4. 

The ratio a/A 2 L of the string tension to the square of the lattice scale 



parameter is hard to measure accurately fllql , and our procedure is subject 
to systematic errors due to our measurement scheme and to the small size 
of our lattice. To minimize these systematic errors, we applied the same 
procedure to the unmodified Wilson action and computed a ratio of ratios: 
the ratio app/K 2 L divided by the ratio aw / 'Af in which app is the positive- 
plaquette string tension, aw is the Wilson string tension, and is the lattice 
scale parameter. 

We performed simulations of SU (2) gauge theory guided by the unmodi- 
fied Wilson action on an 8 4 lattice and measured the ratio aw / 'Af by the same 
procedure that we used to measure the positive-plaquette Wilson action. We 
found a narrow transition region between (5 = 2.125 and (3 = 2.5 in which our 
best fit was (5 e = 0.43 A Anin = 1.42, /5 max = 2.64, and a w /A-l = 30.8 ±4. In 
the interpolation (||), the coefficient (/5 max — /3)/(Anax — Anin) of the Wilson 
string-tension term aw^iPe) runs from 0.42 at /3 = 2.125 to 0.12 at (3 = 2.5, 
and the coefficient — Anin)/(Anax — Anin) of the tree- level perturbative term 
XoihjiPe) runs from 0.58 at (3 = 2.125 to 0.88 at (3 = 2.5. Our fit is shown 
in Figure 5. 

The rejection of plaquettes that are of negative trace is an entirely non- 
perturbative modification of the Wilson action. To check this assertion, we 
ran a positive-plaquette simulation at (3 = 16 and found that the update 
algorithm never tried to form a plaquette of negative trace. Thus for suffi- 
ciently large values of (3, the positive-plaquette Wilson action is practically 
equivalent to the usual Wilson action. The lattice scale parameter A^ is 
defined JTTJ] as the limit 



A L ^ hmiexp (- -LJ) (P g 2 )- M2 ^. (10) 
The value of A^ appropriate for the positive-plaquette Wilson action is there- 



fore the same as that for the unmodified Wilson action |TI| . And so our ratio 
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of the positive-plaquette ratio app/A 2 L = 1.61 ± 0.2 to the Wilsonian ratio 
o"vi//A| = 30.8 ± 4 may be identified as the ratio of the positive-plaquette 
string tension opp to the Wilson string tension aw, 

= —= °- 05 - en) 

Hence somewhat more than 90% of the confinement signal measured in Wilso- 
nian simulations on an 8 4 lattice is due to plaquettes of negative trace, which 
are lattice artifacts. 

Because our simulations were carried out on a small lattice, we make no 
statement about any possible differences that might exist in the continuum 
limit between the physics of the standard Wilson action and that of the 
positive-plaquette Wilson action. In fact we imagine that any such differ- 
ences are negligible, although the approach to the continuum limit may be 



smoother for the positive-plaquette action |L2], pflfl . Indeed beyond ft = 16, 
the Wilson action and the positive-plaquette Wilson action are essentially 
equivalent. But confinement has been clearly exhibited only on lattices of 
moderate size (< 100 4 ) at relatively large couplings with actions that are ei- 
ther compact || [|, |17j or have auxiliary fields || |9| . Thus our result provides 
motivation for further studies of confinement at weaker coupling on larger 
lattices with actions that suppress artifacts, such as the positive-plaquette 
Manton action fl8|1 . Such studies might suggest effective actions suitable for 



use at stronger coupling. 
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Figure 1: The U(l) Creutz ratios x(hj) as given by Wilson's action and the 
exact Creutz ratios, x(2, 2) (solid), x(3, 3) (dashes), and x(4, 4) (dots). The 
confinement signal at j3 — 0.75 is striking. 
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Figure 2: The £/(l) Creutz ratios xihj), as given by Wilson's action with 
a small barrier, and the exact Creutz ratios at (3 e = 0.73/3, x(2, 2) (solid), 
x(3,3) (dashes), and x(4, 4) (dots). There is no false confinement signal. 
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Figure 3: The SU(2) Creutz ratios x(hj) as given by Wilson's action 
restricted to positive plaquettes and the tree-level perturbative formula 
Xo{i,3,Pe) f° r the Creutz ratios x(2, 2) (solid), x(2, 3) (dots), x(2,4) (dots), 
and x(3, 3) (dashes) at the effective inverse coupling (5 e = 0.54/3. 



13 



1.0 




1.0 1.5 2.0 2.5 3.0 

a 

Figure 4: The SU{2) Creutz ratios x{hj) as given by Wilson's action re- 
stricted to positive plaquettes and the fit @ to the Creutz ratios x(2, 2) 
(solid), x(2, 3) (dots), x(2,4) (dots), and x(3, 3) (dashes). 
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Figure 5: The SU{2) Creutz ratios x(hj) as given by the unmodified Wilson 
action and the fit @ to the Creutz ratios x(2, 2) (solid), x(2, 3) (dots), x(2, 4) 
(dots), and x(3,3) (dashes) between (3 = 2.125 and (3 = 2.5. 
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